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Abstract 

The presence of a massive scalar field near a Kerr black hole is known to produce 
instabilities associated with bound superradiant modes. In this paper we show that for 
massive fermions, rather than inducing an instability, the bound superradiant modes 
condense and form a Fermi sea which extends well outside the ergosphere. The shape 
of this Fermi sea in phase space and various other properties are analytically computed 
in the semiclassical WKB approximation. The low energy effective theory near the 
black hole is described by ripples in the Fermi surface. Expressions are derived for 
their dispersion relation and the effective force on particles which venture into the sea. 
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1 Introduction 

A massive scalar fleld has a dramatic effect on a rotating black hole. Small fluctuations 
of the scalar can grow exponentially - the so-called "black hole bomb" [ij - and classically 
destabilize the Kerr solution. Recently it has been argued that this phenomenon both 
constrains the possible spectrum of light axions in our universe and predicts certain gaps in 
the mass-angular momentum spectrum of astronomical black holes. 

The unstable modes arise from bound states with energies below the superradiant bound. 
They are localized in a "pocket" in the effective radial potential which is separated by a 
barrier both from inflnity and the horizon. One way to view the instability (in the particle 
picture) is this: thermodynamically the superradiant modes want to form a Bose condensate 
in the pocket. This condensate then stimulates emission from the black hole, which leads to 
a larger condensate and even more emission. This runaway behavior is the black hole bomb. 

It is not known if there are light scalars in the real world that could lead to such insta- 
bilities. We do know however that there are light fermions. In this paper we investigate 
the effects of a light fermion on a rotating black hole in the semiclassical WKB limit, which 
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turn out to be equally dramatic but qualitatively different. The WKB limit is relevant for 
fermions whose Compton wavelength is small compared to the black hole radius. For a 
neutrino near a solar mass black hole the WKB expansion parameter is around 10~^. 

The leading WKB wave equation for a massive scalar and a massive fermion are identical 
and equivalent to the massive geodesic equation. Therefore the fermions also see a pocket 
outside the black hole with bound superradiant modes. However fermions cannot form 
a Bose condensate. Instead, in the standard Unruh vacuum, all the bound superradiant 
modes are filled, forming a Fermi sea which extends well outside the ergosphere. The low 
energy effective theory is described by ripples on the Fermi surface which can carry energies 
much below the fermion mass. In the WKB approximation, many properties of this Fermi 
sea turn out to be analytically calculable. We compute the number density of fermions, the 
fluid continuity equation and the force on a test particle which couples to the sea fermions. 

In principle, it is conceivable that our results are relevant for understanding the properties 
of astronomically observable black holes, and the Kerr-Fermi sea a component of dark matter. 
In this paper we consider only the simplest case of a free fermion on a fixed geometry, and 
a number of issues must be addressed before our results can be applied to the real world. 
Among these are the effects of fermion interactions, gravitational backreaction and the rate 
at which the quantum state around a black hole formed by gravitational collapse approaches 
the Unruh vacuum. 

In a completely different direction, our results may also be useful in the search for holo- 
graphic duals for rotating condensates for example of the kind created in the laboratory in 
i- 

Recently it has been argued that the near-horizon deep-infrared region of extreme Kerr 
is dual to a two-dimensional conformal field theory (CFT) Q. The Fermi sea described in 
this paper is localized a finite distance outside the horizon, and therefore is largely relevant 
at energies scales higher than those considered in [Q. Nevertheless the sea does touch the 
horizon and so may be relevant to the near-horizon CFT. Fermi seas of this type - localized 
outside the black hole - have been studied in the context of holographic superconductivity 
in 

This paper is organized as follows. In section 2 we review the definition of the Unruh 
vacuum and describe the semiclassical limit. In sections 3 and 4 we find the shape of the 
Fermi surface in phase space, compute the number density and other feature of the ground 
state, and derive the equation governing small fluctuations. In section 5 we compute the force 
on a test particle which couples to the fermions along its worldline. Although we anticipate 
that our results hold qualitatively for general Kerr, for the sake of brevity we restrict to the 
maximally rotating case. The non-extremal case is considered briefly in appendix A and 
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details of the WKB approximation are provided in appendix B. 



2 Semiclassical limit of the Kerr-Unruh vacuum 

A Schwarzschild black hole has two well-known time-symmetric vacua: the Hartle-Hawking 
vacuum, which has a thermal gas of particles and a smooth stress tensor at the horizon, and 
the Boulware vacuum which has no particles (according to static observers) and a singular 
stress tensor at the horizon. In addition there is a third time-asymmetric but stationary 
vacuum, the Unruh vacuum, which approximates the late-time asymptotic state around a 
black hole formed by gravitational collapse. It is characterized by no "in" particles from past 
radial infinity and a thermal flux of "up" particles from the past horizon H . These initial 
conditions result in outgoing Hawking flux at future infinity and a smooth stress tensor at 
. The Unruh vacuum describes a black hole emitting blackbody radiation into empty 
space. 

For the Kerr black hole, analogs of the Hartle-Hawking and Boulware vacua are known 
generically not to exist, see [^, |TD|, 0, |T2|, [T^, |T3|, [T^ for discussion. The basic problem is 



that one cannot equilibrate modes with respect to the rotational potential outside the speed- 
of-light surface. The Unruh vacuum on the other hand does exist. For the sake of brevity we 
will specifically consider only the case of extreme Kerr except for appendix A where a few 
details of the general case are given. In the extreme case the Unruh vacuum is characterized 
by no incoming particles from past infinity and a "thermal" fiux at zero temperature with 
nonzero angular potential VLh = • This populates with unit probability all up 

modes with energy below the superradiant bound: E < itlVLh where m = is the angular 
momentum of the mode. More explicitly, in terms of the expansion of the field operator for 
a Dirac fermion, 

(2.1) 

where i labels the angular eigenfunction, the Unruh vacuum is defined by 

a^"|0) = 6*"|0) = a^^'lO) = ^"^'lO) = (2.2) 

for all a, h. 

The massive Dirac equation cannot in general be solved analytically for Kerr so, although 
we know in principle which modes are occupied in the Unruh vacuum, it is difficult to 
determine quantities such as for example the stress tensor. In this paper we are interested 
in the leading semiclassical or WKB approximation, in which case it turns out the leading 
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behavior of many quantities of interest is analytically calculable. The semiclassical expansion 
parameter is the ratio ^ of the Compton wavelength of the particle to the Schwarzschild 
radius. For a neutrino near a solar mass black hole this is -Arj ~ 10"''. 

The semiclassical approximation begins with the phase space of a massive particle in the 
Kerr geometry which is labeled by the 6 coordinates {r,9,(j),Pr,pe,P(f>) ■ This phase space 
is constrained by the required existence of a real four vector pi, obeying p^Pu = —fJ''^, with 
Pt = —E, and lying in the future light cone.Q We then define single particle states as cells of 
volume in this phase space. A quantum state on the entire spacetime is specified by giving 
the occupation numbers of each of these states.^ For fermions it can only be zero or one. 
For Dirac fermions, as shown in appendix B, we need four copies of this space corresponding 
to spin up and spin down particles and anti-particles. 

Cells in phase space are mapped to orbits by the action of the Hamiltonian, which in 
turn lift to pieces of WKB limits of solutions of the massive Dirac equation. The cell is then 
filled if the corresponding semiclassical field mode is. Since the Unruh vacuum is stationary, 
all cells on the same orbit are filled with the same probability. Clearly a necessary condition 
for an orbit to be filled is 

E < mnff. (2.3) 

In general, the WKB limit of a given solution has real "under-the-barrier" regions as 
well as oscillatory regions. More than one oscillatory region corresponds to more than one 
classical orbit separated by a potential barrier. Moreover, it may happen that one or more 
of the oscillatory regions corresponds to an "anti-orbit" with momentum vector p^ which lies 
in the past rather than the future light cone. Indeed this occurs in flat space. In the case of 
fermions these anti-orbits correspond to the negative-energy particles filling the Dirac sea. 
In the case of Kerr we shall encounter anti-orbits with positive energy inside the ergosphere. 

Orbits are characterized by the three conserved quantities E, m and the Carter constant 
Q (defined in the next section) which are functions on phase space. Let us now consider the 
various possibilities in turn. 

2.1 E>i2 

Orbits with E > fi have energy greater than their rest mass and are therefore not bound. 
They all reach infinity in the far past and/or future. Those which do cross the horizon must 

^This constraint can be expressed as > —gt^m. 

^The notion of a particle depends on the choice of a coordinate system. In this paper we will use Boyer- 
Lindquist coordinates. Our operational definition of a particle in a given orbit is an occupied mode of the 
associated WKB wavefunction. This may become unphysical inside the ergosphere where static observers 
in Boyer-Lindquist coordinates are moving above lightspeed. In section ^, we will compute something more 
physical: the force on a probe which couples to the fermions comprising the Fermi sea. 
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have 

P,xti < 0, (2.4) 

where xh = dt + ^nd^ is the null generator of the horizon, as both p and xh are future- 
directed at the horizon. The condition ( p.4|) may be rewritten d& E > mVlH-, so we conclude 
that such unbound orbits are unfilled. 

\i E > ^ and E < rnQn, the effective potential governing the WKB wavefunction has 
a barrier separating the near horizon region from the far region, under which the wave 
function must tunnel. In the semiclassical limit, this barrier is infinite, and the eigenmodes 
on the near and far side of the potential do not mix. Those on the far side begin at spatial 
infinity and return there after reflecting off the barrier. Such modes are not populated in 
the Unruh vacuum. Those on the near side begin and end at the horizon. However if they 
have E < mVlH and cross the horizon, the momentum vector must lie in the past light cone 
and they correspond to anti-orbits. These modes are filled in the Unruh vacuum in the same 
sense that the negative energy modes in flat space are filled in the usual Dirac sea. 

Away from the strict semiclassical limit, the WKB wave functions on the near and far 
regions have exponentially suppressed tunneling through the potential barrier. This is the 



well-known Unruh-Starobinsky radiation of fermions [|T7 



2.2 E < tuQe and E < fi 

These orbits cannot reach infinity because their energy is less than their rest mass. They 



also cannot cross the horizon as they do not satisfy (|2.4|) . They are bound orbits which 
remain a finite distance outside the black hole. One expects that they are filled since they 
are near the horizon which behaves like a reservoir for modes with E < ttiVLh- We will see 
this is indeed the case. 

To illustrate the point, let's take a specific orbit on the equator with say E = ^ and 
m = ^. In the WKB limit the radial equation for the wave function becomes a Schrodinger 
problem 

+ V)^P = 0, (2.5) 

with a potential 

where s = log (r/M — 1), illustrated in figure |1|. 

This potential is of course the same as the one that appears in the geodesic equation for 
the corresponding orbits. The horizon is at r = M, s = — oo. V goes to — oo at the horizon 
and +00 far from the black hole. In between there is a local minimum at Smin = 0.93, and a 
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Figure 1: Potential at E = ^fi, m = 2Mfi in units of iJ?M'^. The horizon is at s = — oo and 
the boundary of the ergosphere at s = 0. 

maximum at Smax = —1-97. The zeros are at {si, 82,83) = (—2.61,-0.27,1.55). The WKB 
solution is 

z^(s)~e-/^'^^'v^, (2.7) 

where we take the branch of the square root which is positive at s ^ 00 where the wave 
function vanishes. At s ^ —00 ip ~ e^*^^ "-1-^1+1211^14) ^ Since the wave function emanates 
from H~ rather than radial infinity and has E < mVlH the mode is filled in the Unruh 
vacuum. 

The zeros of V (51,52,^3) divide the real and oscillatory regions. The two oscillatory 
regions have —00 < s < Si and S2 < s < S3. In the first region g^ci,E < —gt^m, so we have an 
anti-orbit. The filled mode does not correspond to a real particle in the sense that a geodesic 
observer very near the horizon will not see any particles. Note that even though E, m and Q 
are the same everywhere, it need not be the case that the tangent vectors in both oscillatory 
regions are both in the future light cone. The second region indeed has g^p^^E > —g^m and 
corresponds to a bound orbit in the equatorial plane (since we took Q = 0) around the black 
hole. Since this orbit arises as the WKB limit of a filled mode, all the fermion states in this 
orbit are filled. 



We note that the so-called Penrose orbits, which have E < 0, can readily be shown |jT6 
to always have E > itlVLh- Hence filled orbits always have positive energy. 

To summarize, the filled portion of the Fermi sea is the region of phase space obeying 

E < min(/i, mfiji/). (2.8) 

These are the bound orbits with energies below the superradiant bound. We will refer to 
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the E = mVtH surface as inner and the E = jj, surface as outer because the typical radii of 
orbits grows with angular momentum. 



3 The inner Fermi surface: mVtH < 

The Fermi sea and surface each have two portions divided along the 5-dimensional hyper- 
surface m = ^/^H in phase space. In this section we analyze the inner surface m < ^/^h- 
The outer surface with m > fi/Qn is analyzed in section ^. 

3.1 The ground state 

The extreme Kerr metric is 

ds^ = -4 (dt - Msin^ ed(f)? + ((r^ + M^)d(f) - Mdt? + ^dr^ + P^dO^ (3.1) 

A = (r - M)2 , p2 ^ + cos^ 6 , 

where we take Planck units G = h = c = 1. The horizon is at r = M. The angular velocity 
of the horizon is 

= (3.2) 

The 6-dimensional phase space for a particle in the Kerr geometry has coordinates (r, 6, 0, pr, pe, fn) 
where we denote m = p^. The Fermi surface is a 5-dimensional hypersurface in the phase 
space, given by 

E{r,9,(l),pr,pe,m) = rnQn (3.3) 
where, for particle mass fi, E = —pt is the solution of the quadratic equation 

P^P, = -/^' (3.4) 

such that is in the future light cone. It is sometimes convenient to eliminate pe in favor 
of the conserved quantity WE 



_ (m - ME sin^ ^)^ 999./ .9 , . 

Q = pI + ^ 5 ^ + MV cos^ ^ - (m-ME)2 (3.5) 

sin 9 

known as the Carter constant. In terms of Q the general form of potential appearing in the 
radial equation Ap^ + \^ = 0, or equivalently in the WKB Schrodinger equation ( p.5|) is 



y = -^ {E{r^ + M^) - mMf + /iV^ + (m - MEf + Q. (3.6) 
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The insertion of the surface condition (|3.3|) into (|3.4|) yields 



Api = xm'-pi-fi'p' (3.7) 

r2 + 2Mr 2 2 2 n 
m — jj, r — Q 



where 

X = -^^\d, + nj,d^\' (3.8) 
A sm 6 

= - H \ - csc^ 9 sm^ 9 

4 2M 4M2 4 

Therefore for the ground state the Fermi surface is at the real values of 



; - /i2p2 _ p2 
Pr = ±y = ±CTo (3.9) 

where the ± denote the two branches of the hyperbola. Note that near the north pole, 
sin^ —>■ 0, and x ~^ ~oo so Pr can never be real. Hence the Fermi surface does not extend 
all the way to the poles. 

3.1.1 Number density 

Its interesting to compute the fermion density. The measure on phase space is, for h = 1, 

(ir Ad9 Ad(p A dpr A dpg A dp^. (3.10) 

The total number of fermions present is just the volume of the Fermi sea computed with this 
measure. The density of fermions as a function of (r, 9,(f)) is 



N{r,9,(j)) = - — —-= / dprAdpeAdp^ (3.11) 



sea 
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(27r)3v^ 



2Mfi /" xm^— /i^p^ 



dm / dpgO-Q 



{2TiYVhK^ 3 3yx' 



where h is the induced metric on a constant time slice and 



r- ^ / (r2 + M2)2 - M2A sin2 ^ 
Vh = psm9d^— . (3.12) 

For radii r of order M, p ~ M, A ~ M'^, \fh ~ NP' and x ~ 1- It then follows that 
N fi^. This is just the statement that in the vicinity of the black hole there is of order 
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one fermion per Compton cube. At large r, we have p ~ r, A ~ ~ r and x ~ JW^- 

In this hmit the range of the m integral goes to zero and the leading and subleading terms 
in cancel. The leading nonvanishing term gives ~ There is a log divergence 

at large r in the total number of fermions. 

3.1.2 Angular momentum density 

It is easy to compute the angular momentum density L by just inserting a factor of m in 



the Fermi sea integral (|3.11|) . One finds 



^(r, 0) = "".f rj-r {^M\ - (3.13) 

Again there is a large r cancellation, and the leading behavior is L ^ jj, ^M{^f. The angular 
momentum of the black hole itself is J = M^. If we impose a large radius cutoff Tc, the ratio 
of the Fermi sea angular momentum to the back hole angular momentum isQ 

^-/M^ln^. (3.14) 

In Planck units, the mass of the sun is 10^^, while a typical neutrino mass is 10^^^. 
Hence for the neutrino sea around a solar mass black hole, for L to be of order J we need 
^ ~ 10^°"*. For r of order M, j ~ 10"^^. In this example the semiclassical expansion 
parameter jjj ~ 10^^ is very small. The computation of the energy density is a bit more 
involved but one expects on dimensional grounds E ~ which is also small. Hence 

the backreaction of the Fermi sea on the geometry is suppressed. 

3.1.3 (^'^) 

In the WKB approximation (discussed in appendix B), the current for a single particle mode 
ip± labeled by is 

^7/^^ = ±-'^P'^|Jp^' + 0{h) , (3.15) 
/i 

where the sign is positive for particles and negative for antiparticles. Because the sea has an 
equal number of particles and antiparticles, the total current vanishes 

^ijr^M^) = . (3.16) 



■^The total angular momentum also has a logarithmic divergence at the horizon. This divergence does 
not appear away from extremality, and is similar in form to well-studied divergences of the stress-tensor 
near the horizon of an extremal Reissner-Nordstrom black hole. In the latter case the divergence appears to 
be rather benign and eliminated by gravitational backreaction - see jlH and references therein. A similar 
outcome may pertain to extreme Kerr. 
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In the semiclassical approximation, has real and positive contributions from both 

particles and antiparticles. The expectation value is given by summing ( |3.15| ) over the Fermi 
sea, 



4/i 



(27r)3 
Svr/ip 



dmdpodpr 



9 J sea 

2 rpi-I^H 



1 



(3.17) 



dm {mflH - Erain) 



where the factor of 4 comes from particles and antiparticles of each spin, and 

1 



Emin = - { 2MVm + p^/ A(p2mV sin^ 6) + Aa/i^ 



with 



9t^ 



(3.18) 



(3.19) 



(r^ + M^)^ - AM^ sin^ 6 
is the minimum allowed energy of an orbit given r, 6, 0, m. Doing the final integral we find 



(27r)2v/=^ 



sin log I — + + 



S'tfli ^ mp VA /m^ a sin 6' 

p^ apsin^^ y p2 p2 



(3.20) 



' m 



p 



p- 



a sin 9 



3.2 Small fluctuations 

In this section we derive the first-order differential equation governing the propagation of 
small disturbances of the Fermi sea. Denoting the top of the sea by a = Pr{r,6,(f),p^,m), 
the general equation for the evolution of the Fermi sea is 



dta = [pr, H] + dp^a[p^, H] + Sp.afpe, H] - dr(j[r, H] - d^a[(l), H] - dea[9, H] 
where p^ = m and H = E. The second term vanishes. Now we expand 

a = ±ao + a^. 
( p.21| ) then becomes the fluid continuity equation 



which is related to conservation of \I'7'^\I^. Here 



dH 

dpi 



(3.21) 



(3.22) 



(3.23) 



(3.24) 
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A^r2 + 2Mr-^(/xV2 + Q) 

r(r2 + Mr + 2M^) + M^{r - M) cos^ 6 

2M(r - M cos^ csc^ 9 
r{r^ + Mr + 2M2) + M^{r - M) cos^ 9 

2M(r - M)^Q- cos2 ^[/i2M2 + m2(-3^ - 



m r 



r2 + Mr + 2M2) + M2(r - M) cos2 ^] 



(3.25) 
(3.26) 

(3.27) 



4 The outer Fermi surface: m^H > 

In section 3 we found the inner Fermi surface for < given by the condition E = mQu- 
Here we describe the case mflH > fi in which case the Fermi surface is at ii^ = yU where the 
orbits become unbound. In phase space the equation E = fi becomes 



Pr 
C2 
Cl 

Co 



±0-0 = ±— ycam^ 
M^ - A csc^ 9 
ArM^ 

2Mr(r^ + M^) 



ficim + yu2co — pjA 



(4.1) 



Outside the ergosphere C2 < and the resulting contribution to the number density is 

4 



N 



i'fj.p+ p^/ {c2m^-ficim+iJ.^co)/A 

(O \Z fh I "^"^ '^P'''^ 
[Znj-^Vri J2M^l Jo 



2co(p+-2M)+Ci(4M2-j9^ 



P+ 



2(27r)3v^A 
sin 9 (^-2M'^r sin 9 + V2Mr{r - M)p 



^-fm^-pl) 



(4.2) 
(4.3) 

(4.4) 



(r2 - 2Mr + M2 cos2 9) 
At radius of order M, the number density is also of order /x^ as in (|3.11| ) so N N. However 



at large radius ( [4.2| ) falls off more slowly than ( |3.11 

N 



,M. 



^ — ^: 

r 



as r 



oo 



(4.5) 



The angular momentum density goes as L ~ /i^M(^), while the energy density goes as 
E ~ fJ^^i^)^^'^- For r ~ M, the number, angular momentum and energy densities are 
comparable to those of the inner sea. However they fall off more slowly at large r. 

The computation of (^\I') is identical to that of the inner Fermi sea, except that in the 



final answer ( |3.2CI| ) the limits of integration are m = ii/VLh to m = ^ip+. 
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5 Sailing on the Fermi sea 



In this section we consider the effect of the Fermi sea on the trajectory of a particle with 
worldline couphng 

A / dr^'^H (5.1) 



to first order in the couphng A which we take to be smah. Here r is the proper time along 
the particle trajectory X'^(r), which implies g^i.drX'^drX'^ = — 1. There are two effects to 
consider: the change 6Ep in the energy of the Fermi sea and the change 6Ep in the energy 
of the test particle. Both of these result in an effective force on the test particle, and a 
modification of the geodesic equation if the particle is in free motion. We will assume in this 
section that the particle motion is slow enough so that effective description of the Fermi sea 
in terms a given by ( p.23| ) and the semiclassical description are valid. 

The change in the particle energy is read directly off of the interaction ( p.l| ). It is 
nonvanishing at leading order in A because has a nontrivial vev, as computed for the 
inner sea in section 3. One finds simply 

dr - 

6Ep = X—{^^) (5.2) 
To compute SEp, note that the presence of the particle modifies the Dirac equation to 

- (/i + = (5.3) 

with 

6fi = f dT6\x^ - X^{t) = -^'^6\x' - X\t)) (5.4) 

V -fi' J y-9 dt 

This in turn modifies the equation for the Fermi surface, which at linear level can be expressed 
as a source term for a. The modification is found simply by replacing /i by /i + 5/1 in (|3.9|) . 
This gives for the inner sea 

= T , '^^^\x' - X\t)) (5.5) 

V (-^) A(X"^^ - /"V - Pe) 

The contribution to the energy of such a "dimple" from every m and pe is 

Am/ip^ dr 

MV(-^7)A(xm2-pV-P^) 

Now we integrate over and m and multiply by 4 (for particles and antiparticles of each 
spin) to get the total energy. This gives, where x is positive, 

^ 2A.V(r-2Mcot-.)d. 
(27r)2xv/F^ dt 
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The gradient of this gives a force which vanishes at infinity and diverges on the horizon. At 

3 3 

large r, SEp goes as ^ and the force as The sign depends on the sign of A. 
There is also a dimple in the outer sea. Outside the ergosphere the energy is 

4AyU^p^ (sin 9V2Mrp - 2M(r - M cos t/j ] 
^ (27r)2^/=^(r2-2Mr + M2cos2 0) dt ^ ' ^ 

At large r, 6Ep goes as /i'^r~i, and the force goes as /i^r~t. 

We note that the coupling ( |5.1| ) has the special property that, at first order, the presence 
of a particle only produces a dimple in the Fermi surface which follows the particle trajectory, 
but not a propagating wave. For more general couplings - such as J ^jf^'^dX^ - we do expect 
waves to be produced and propagate according to the fluid continuity equation ( |3.23| ). 
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A Non-extremal Kerr 

We now consider a black hole at temperature Th, with rotational parameter a = J/M, 
horizon at r = r+, and 

A = r^ + a^-2Mr, = + cos^ 6 , = , " o • (A.l) 

The inner Fermi surface, defined by setting E = uiVLh in the geodesic equation, is 



where 



Pr = ±\i^^^^^-^^^^ (A.2) 



X = -r^Adt + nud^? (A.3) 
A sm 6 



[^//(r^ + a^) - a]2 (l-afi/^sin 



2 Q\2 



A sin^ e 

At finite temperature the Fermi surface is smoothed out, but modes with 



(A.4) 



Q<E<mVlH, \E~mVlH\^TH (A.5) 
are filled with unit probability. 
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B WKB limit of the Kerr-Dirac equation 



The WKB limit of the Dirac equation in a general curved spacetime is considered, for 
example, in pO| , Expand the fermion wavefunction as a series in h , 



^ = exp(-^A/n) ^{-ihy 



ttrAX 



(B.l) 



n=0 



Plugging this into the Dirac equation, we have the leading order equation 

- /i)ao = . 

The existence of a nontrivial solution for requires 



fB.2l 



(B.3) 



Therefore = —A^^ is a congruence of timelike geodesies pO| , |21| . Focusing on a single 
worldline of the congruence, for modes moving forward in time we can adjust the local 
tetrad so that 

Pa = ej:p^ = (-/i, 0,0,0). (B.4) 
Then the leading order positive energy solution is 

= gi/<ixMp, (^^^ 0, 0)^ 



along the worldline, where 

To leading order in h, the fermion bilinears are 



I |2 I I |2 



■A* 
J+ 



P' 



where we have used the Gordon decomposition to compute the current. 
Similarly, the negative frequency solution is 

^_ = e-hf<i-'p. (0,0,/?i,/52) 



(B.5) 
(B.6) 

(B.7) 
(B.8) 

(B.9) 



(B.IO) 
(B.ll) 



*Our gamma matrices are 7* = 



cr' 

-a' 



7*^ = ( Q ) , where ct* are Pauli matrices. 
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with the bihnears 



= + (B.12) 

So far, we have shown that to leading order the Dirac equation decouples into four 
components, particles and antiparticles of spin up and spin down, where each component 
travels along an identical particle geodesic. We now specialize to the extreme Kerr metric. 
There are three conserved quantities, the energy E, the angular momentum m, and Carter's 
constant Q. The Dirac field can be written as a mode sum 

dEdmdQ {aip^ + b^tlj-) (B.14) 

where the WKB wavefunctions ip± and a, b depend on E, m, Q, s. Comparing to (|2.1| ), we 
see that the Unruh vacuum has 

a\0) = b\0) = (in modes) (B.15) 
a|0) = b\0) = (up modes, E > mVtH) 
a^|0) = 6^10) = (up modes, E < mfi/^) 



Then we have 



= I dEdmdQ {j1 - j^) (B.16) 

J sea 

[^'^)=[ dEdmdQ^^{iX+ 1) (B-17) 

J sea P 
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